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■ Abstract. Three case studies are presented of slepton pair production followed by 
two-body or quasi- two- body decays at a muon collider. Precision mass measurements 
are possible using a variety of kinematic fitting methods. Standard Model and super- 

j-^. symmetric backgrounds are easily controlled by kinematic cuts. In all three cases it 

£S) . appears that detector resolutions, not backgrounds or statistics, will dominate the final 

error bars. Polarized beams are not necessary to control SM backgrounds. However, 
£^ . without polarization it may be difficult in some cases to disentangle Ir from li signals. 
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^ Introduction 

A muon collider is in principle an excellent machine for precision studies of weak 
scale supersymmetry. Depending on y/s and the SUSY mass spectrum, it may 
be possible to observe pair production of a half-dozen or more distinct sparticles. 
For R parity preserving SUSY, sparticle pair production is kinematically under- 
constrained, due to the pair of unmeasured LSP's. However in many cases each 
sparticle in the pair has a significant branching fraction for what is essentially a 
two-body decay: 

sparticle — > LSP + particle , (1) 

where "particle" refers to a fully reconstructible Standard Model particle (e, W, 
Z, and possibly h , t), while the LSP is assumed to be the lightest neutralino Xi- 

In these cases there are a variety of kinematic fitting methods for extracting 
sparticle masses. In this talk I will report on two such methods applied to smuon, 
selectron, and sneutrino production at a muon collider. Chargino production is not 
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discussed, since for light fermionic sparticles the best method for a precision mass 
measurement is a threshold scan [1]. An interesting challenge for future investiga- 
tion is the production of staus, stops, and the heavier chargino and neutralinos. 

Sparticle production at a muon collider is similar in many respects to sparticle 
production at an e + e~ machine. For the present analysis the most important differ- 
ences are that the muon collider has (i) much higher energy reach, (ii) significantly 
lower advertised luminosity at comparable energies, (iii) little or no polarization 
available without taking a significant hit in luminosity, and (iv) large detector 
backgrounds from muon decays. These detector backgrounds are generally soft, 
but large fluctuations could cause problems for precision SUSY measurements. 
They will also impact on isolation cuts, determinations of missing Et, and detector 
resolutions generally. These problems will be left to future study. 

At a muon collider smuon pairs arise from both s and t channel production; the 
s channel production is through a virtual photon or Z, while the t channel diagram 
involves the exchange of a neutralino. The s and t channel contributions interfere 
destructively, but this effect will not be important for the examples considered 
here, where the t channel production is dominant. Selectron production proceeds 
only through the s channel, and is thus suppressed in the examples. Muon sneu- 
trino production proceeds only through the t channel, and is thus competitive with 
smuons. 

In both supergravity (sugra) and gauge mediated models, the Irs are lighter 
than the Il's. Independent of the SUSY model, Irs decay almost 100% via a single 
two-body mode: 

fiR->XiP I e R -> xi e • ( 2 ) 

The branching fractions of the l^s and z> L 's are model dependent. The important 
decay modes for, e.g., p,L are: 

xf ^ 
Kinematics 

The basic kinematics can be understood by considering pair production of /2r: 

Vr(P2) -> X°i(P5)KPe) ■ (4) 

Each event consists of an acoplanar dimuon pair plus missing Et- Six measurements 
are made, i.e., the 3-momenta of the two muons. The event is characterized by 13 
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TABLE 1. Sparticle and Higgs spectrum for LHC Point 5, 
which corresponds to minimal sugra parameters mo=100 GeV, 
m 1/2 =300 GeV, A o =0, tan/3=2.1, and sgn(/x)=l. 



Particle Mass (GeV) Particle Mass (GeV) 
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kinematic variables: the four 3-momenta of the final state plus the common LSP 
mass. There are 5 kinematic constraints: one from the assumption that the two 
smuons have the same mass, and the rest from the known initial state 4-momentum. 
This leaves 2 undetermined variables in the event, which we may take as M^, M LSP . 

The kinematic endpoint method arises from the following expression for the en- 
ergy of each muon as measured in the rest frame of its parent smuon: 

E o _ mi - M L 2 SP 

where we are neglecting the muon mass. The maximum and minimum boosts from 
this frame to the lab frame then provides us with two kinematic endpoints E™ ax , 
E™ m , in the muon energy spectrum. A precision measurement of both endpoints 
allows us to extract both and M LSP . Note that this method requires good 
statistics to be useful, and does not take advantage of all the kinematic information 
in the event. 

Another kinematic method, developed by Feng and Finnell for e + e~ studies, 
extracts assuming a precision value of M LSP is already known from other sources 
and thus can be used as an input. This method starts with the relation 

M? = - \p 3 \ 2 - \pi\ 2 - 2\p 3 \\p 4 \cos 9 M . (6) 

For given input value of M LSP , the only unknown on the right hand side is # 34 , the 
angle between the 3- vectors p 3 and p 4 . This angle is then estimated, event by event, 
by a certain function of measured variables. This function has the property that 
the error of the estimate goes to zero in the limit that the two LSP's are back-to- 
back in the lab frame. For \fsj1 ^> ^> M LSP , the mass estimates peak strongly 
around the true value, and precise results are possible even for rather sparse data. 

A third kinematic method, which is currently under investigation, involves adapt- 
ing the likelihood methods developed for extracting the top quark mass from the 
dilepton channel. This method is also well-suited to sparser data sets. 
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Simulations 



Simulations were performed using PYTHIA v6.1 [3] coupled to the ATLFAST 
vl.25 [4] fast detector simulator. Note that the small differences in the sparticle 
spectra produced by PYTHIA and ISAJET [5] make a difference for the analysis 
done here. The ATLFAST defaults were used for lepton isolation and jet recon- 
struction. Smearing was not included, and no attempt was made to include detector 
backgrounds. Thus "precision" here refers only to statistics and to SUSY signal 
versus Standard Model (SM) backgrounds and SUSY backgrounds. 

Sleptons at LHC point 5 

This first study overlaps with the analysis presented by Frank Paige at the Fermi- 
lab workshop [6]. LHC point 5 is a mimimal supergravity reference point described 
in Table 1. For dimuon and dielectron production at t/s=600 GeV, cuts were 
imposed similar to those of [6]: 

• Exactly two isolated e or \i leptons and no jets, 

• E > 10 GeV and \r)\ < 1.3 for each lepton, 




FIGURE 1. Dimuon production after cuts, 20 fbr 1 at V«=600 GeV for LHC point 5. The solid 
line is the total smuon signal. The dashed line is the sum of the Standard Model backgrounds; 
the dotted line is the background from chargino pairs. 



• A0 1>2 < 0.97T, 

• \Pt,i + Vt,2\ > 10 GeV, and 

• missing E T > 20 GeV. 

Note that missing E T signatures are degraded at a muon collider detector, due to 
the 20 degree forward and backward dead cones needed for shielding. This is not 
a crucial point for the present analysis, however. 

The signal acceptance with these cuts is approximately 40%. The cuts are very 
efficient at eliminating backgrounds. The simulations included the six most impor- 
tant SM backgrounds; these are: 

• W + W~ pair production, 

• 7// -> Wu^, 

• Drell-Yan, 

• 111 -> Zfi, 

• 77 — * l + l~ , and 

• ZZ pairs. 




FIGURE 2. Dielectron production after cuts, 20 ftr 1 at V«=600 GeV for LHC point 5. The 
solid line is the total selectron signal. The dashed line is the sum of the Standard Model back- 
grounds; the dotted line is the background from chargino pairs. 



The main SUSY background is from chargino pair production, with both charginos 
decaying leptonically. 

Figures 1 and 2 show the dimuon and dielectron event rates plotted versus muon 
or electron energy, with 5 GeV bins. The SM backgrounds after cuts are rather 
flat and encouraging small, even for the dielectron case. The SUSY background 
is negligible. For 20fb~ 1 of integrated luminosity, Figure 2 also reflects the rather 
poor statistics of selectron production. This is not surprising given that the total 
cross section is only 64 fb. The situation is noticeably better for smuon production, 
where the cross section is 400 fb. 

Figure 3 shows the /x/z— ee flavor subtracted slepton signal, after cuts, broken 
down into the its three components: RR, RL+LR, and LL. The integrated lumi- 
nosity is 100 fb" 1 to enhance the statistics. As discussed in [6], this figure shows a 
rather complicated structure, reflecting the fact that there are eight distinct kine- 
matic endpoints affecting the distribution. These are: 

IrIr : 118 GeV, 9 GeV; 

l R l L : 105 GeV, 11 GeV; 

l L l R : 208 GeV, 40 GeV; 

!J L : 181 GeV, 46 GeV. 
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FIGURE 3. Flavor subtracted slepton signal, 100 fb" 1 at Vs=600 GeV for LHC point 5. The 
solid line is the total smuon + selectron signal. The dashed, dotted, and dot-dashed lines are 
from IrIr, IrIl+IlIr, and IlIl, respectively. 



Comparing with Figures 1,2, it appears that with 20 fb" 1 and a perfect detector, 
one can determine the endpoints at 118, 208, and 181 GeV to an accuracy of one 
bin or better. The other endpoints look very challenging. 

The situation improves if we include the Feng-Finnell estimate for the smuon 
mass. This is shown in Figure 4, plotted with 1 GeV bins. The SM background 
shown is completely negligible. Because of the strong peaking, which actually 
resembles a sharp edge, it is trivial to extract the £lr mass with an accuracy of one 
bin or better. This assumes that the Xi mass is already known to within 1 GeV. 
Similar results are obtained for the e#, with somewhat worse statistics. 



The second study is for the heavy sugra point described in Table 2. The results 
are for dimuon and dielectron production at y/s = 1400 GeV, using the same cuts 
as in the previous example. 

Figure 5 shows the flavor subtracted slepton signal corresponding to 1000 fb -1 
of integrated luminosity. Comparing with Figure 3, one notes several differences. 
In the present case the signal is completely dominated by RR production. This is 
because the branching fraction for /i^ or decay to muon or electron plus Xi is 
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FIGURE 4. Dimuon production after cuts, 20 fb _1 at ^/s=600 GeV for LHC point 5. The solid 
line is the total smuon signal, plotted versus the Feng-Finnell estimate for the smuon mass. The 
dashed line is the sum of the Standard Model backgrounds. 



TABLE 2. Sparticlc and Higgs spectrum for the 
heavy sugra point, which corresponds to minimal sugra 
parameters to =500 GeV, mi/ 2 —350 GeV, ^4o=0, 
tan/?=2, and sgn(yu)=-l. 



Particle Mass (GeV) Particle Mass (GeV) 
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only 16%. At this heavy sugra point, the \xl decays predominantly to either \i u n 
or X2 1 1 - Subsequent decays in these modes are unlikely to pass the cuts. 

Since RR production now dominates, there are effectively only two kinematic 
endpoints: 539 GeV and 106 GeV. Note that the lower endpoint is now sufficiently 
large not to be distorted or hidden by the cuts. Both edges are very sharp in Figure 
5. The SM backgrounds after cuts are negligible. Thus with a perfect detector one 
could extract the masses of both and \i with an accuracy better than 5 GeV. 
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FIGURE 5. Flavor subtracted slepton signal, 1000 fb _1 at - v /i=1400 GeV for the heavy sugra 
point. The solid line is the total smuon + selectron signal. The dashed, dotted, and dot-dashed 
lines are from IrIr, IrIl+IlIr, and IlIl, respectively. 




Figure 6 shows the Feng-Finnell plot for the heavy sugra point. The SM back- 
grounds shown are negligible. Again we see strong edgelike peaking around the 
actual jljt mass of 519 GeV. It is clearly possible to extract the mass with an ac- 
curacy of one bin or better. This assumes that the Xi mass is already known to 
within 1 GeV. Similar results are obtained for the en, but with poor statistics. 

Sneutrino pair production 

Muon sneutrino pair production fits our kinematic scenario, provided that the 
sneutrino has a substantial branching fraction to xf A 4 - The chargino will decay 
predominantly to Xi pl us jets. Thus the signature is an acoplanar dimuon pair plus 
missing plus jets. Note that, in the presence of the R parity violating coupling 
LLE, s-channel resonant production of single sneutrinos may also be possible at a 
muon collider [7]. 

Here we have studied the sugra point described in Table 3, for production at 
y/s = 800 GeV. The branching fraction for the muon into xf A 4 is 56%, while 
the branching fraction for xf Xi pl us jets is 65%. We will employ the same 
cuts as previously, except that we now require two or more reconstructed jets (cone 
radius R = 0.4). 
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FIGURE 6. Dimuon production after cuts, 100 fb^ 1 at y/s—1400 GeV for the heavy sugra 
point. The solid line is the total smuon signal, plotted versus the Feng-Finnell estimate for the 
smuon mass. The dashed line is the sum of the Standard Model backgrounds. 



TABLE 3. Sparticlc and Higgs spectrum for the 
third sugra point, which corresponds to minimal sugra 
parameters m n —225 GeV, mi/2=200 GeV, ^4o=0, 
tan/3=2, and sgn( / u)=l. 



Particle Mass (GeV) Particle Mass (GeV) 
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Figure 7 shows the Feng-Finnell mass estimate after cuts plotted in 1 GeV bins. 
Shown is the total signal from all SUSY production mechanisms. The SM back- 
ground after cuts is negligible. The signal acceptance for muon sneutrino pairs 
after cuts is about 4%. Thus, despite a rather large cross section (over 500 fb) the 
plot has rather poor statistics. Nevertheless we again see strong edgelike peaking 
at the true vl mass of 262 GeV. 

It is interesting to note that a previous study of sneutrino pair production at 




FIGURE 7. Dimuons plus two or more jets, after cuts, 20 fb -1 at -^=800 GeV for the second 
sugra point. Shown is the total SUSY signal, plotted versus the Feng-Finnell estimate for the 
sneutrino mass. 



e + e~ colliders [8] relied on the trilepton plus missing E T plus jets channel to kill 
SM backgrounds. For our study point this does not appear to be necessary. Fur- 
thermore, our complementary dilepton channel has five times the rate, before cuts, 
as the trilepton channel. 

Conclusions 

A variety of precision sparticle mass measurements are possible at a muon collider 
using kinematic methods such as those discussed here. Polarized beams are not 
necessary to control SM backgrounds. However, without polarization it may be 
difficult in some cases to disentangle Ir from II signals. 

It seems likely that in most cases detector resolutions, not backgrounds or statis- 
tics, will dominate the final error bars. Thus it will be crucial to perform simulations 
with a realistic mock-up of a muon collider detector. 

Adequate statistics for the type of analysis presented here correspond to inte- 
grated luminosities of at least 20 fb _1 for y/s in the range 500 to 800 GeV. For 
heavy sparticles and y/s > 1 TeV, the minimum useful integrated luminosity is 
about 100 fb _1 . 

This research was supported by the Fermi National Accelerator Laboratory, 
which is operated by Universities Research Association, Inc., under contract no. 
DOE-AC02-76CHO3000. 
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